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Abstract 

We have studied dynamical behaviour of the infinite-range Ising 
spin glass model with p-spin interaction above and below the transi- 
tion into the non-ergodic phase. The transition is continuous at suffi- 
ciently high external magnetic field. The dynamic critical exponent of 
the power-law decay of the autocorrelation function at the transition 
point is shown to decrease smoothly to zero as the field approaches the 
"tricritical" point from above; at lower fields the transition is discon- 
tinuous. The slow cooling approach is used to study the nonergodic 
behavior below the transition at zero external field. It is shown that 
the anomalous response function A(t, t') contains 5-function as well as 
regular contributions at any temperature below the phase transition. 
No evidence of the second phase transition (known to exist within the 
static replica solution of the same model) is found. At lower enough 
temperatures the slow cooling solution approaches the one known for 
the standard SK model. 



* Submitted to ZhETF 
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1 Introduction 



Free energy surface of the spin glasses in the low temperature phase has a 
very complicated structure and consists of the exponentially large number of 
valleys with infinite barriers between them (at least in the mean field approx- 
imation). This leads to the unique dynamical properties. On the arbitrary 
long (but finite) time scales such systems occupy only one valley and physical 
quantities depending only on spin variables at a single time differ from their 
static values derived from the Gibbs distribution. Usually such a behaviour 
is referred to as a "nonergodic" one. The term "nonergodicity" means here 
that the values of measurable physical quantities (magnetization, suscepti- 
bility, etc.) cannot be considered just as some functions of the temperature 
and magnetic field at the measurement point (T, H); rather, they are junc- 
tionals of the trajectory on the (T, H) plane which lead to the final state one 
is measuring. The quantitative theoretical approach suitable for the study 
of such an nonergodic behaviour in a classical Sherrington-Kirkpatrick (SK) 
spin glass [1] and called slow cooling theory was invented in the end-80's 
by Ioffe et al [2, 3] (see also [4]) and then developed (in a slightly different 
version) in [5]. 

The simplest and most known example of the nonergodic behavior is the 
difference between so-called zero-field-cooled (ZFC) and field-cooled (FC) 
susceptibilities (xzfc and xfc below), which is also well-described by the 
static replica-symmetry-breaking approach by Parisi [6, 7]. However, it was 
shown in [2, 4] that the values of the FC susceptibilities obtained within slow 
cooling approach differ from the Parisi theory results; moreover, the same 
applies even to the values of the internal energy (which might be considered as 
a quite robust quantity). Formally one can understand the origin of difference 
between the results of slow cooling and equilibrium theories as being due 
to non-commutativity of two limiting procedures: the thermodynamic limit 
N — y oo and stationary-state limit t a — > oo (here t a is the time system spend 
in the glassy phase, i.e. ageing time). Within the slow cooling theory it is 
presumed that the limit N — > oo is taken first, whereas in the equilibrium 
theory the limit t a — > oo is presumed to be done before it. As a result, the 
transitions between different valleys of the spin glass phase space which are 
separated by "infinite" (in the limit N — > oo) barriers are strictly prohibited 
within the slow cooling approach, so the contributions of different valleys do 
not follow Gibbs distribution - contrary to the case of the equilibrium theory. 
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As a result, there are two sources of the differences between the values of the 
same physical quantity (e.g. energy) calculated within these two approaches: 
i) the metastable state (valley) within which the system is typically stuck 
after slow cooling is non-optimal, i.e. has higher free energy than the ground 
state; ii) in the equilibrium there is an additional contribution from the sum 
over different valley ( a similar quantity is sometimes called configurational 
entropy or complexity [8]). It is not quite clear to us which of the above 
two sources contributes more; however, there is a clear parallel between the 
temperature dependence of the complexity of the SK model [9] 

S con (T)/N<x(T c -T) 6 . (I) 

and the relative difference between slow cooling and equilibrium free energies 
which behaves in the same way near T c (cf. [2, 3]). Therefore it looks natural 
to expect that either the second above-mentioned reason (item ii)) is the 
decisive one, or both of them are of the same importance. 

It is seen from I that in the vicinity of the spin glass transition (which is 
the only analytically tractable region) the disagreement between the results 
of two theories for physical quantities of the SK model is very weak; indeed, 
relative difference is found to be of the order (T c — T) 3 in the \fc and of the 
order (T c — T) 5 in the internal energy, and did not receive much attention. 
The same statement applies to other spin-glass models that are characterized 
by continuous Parisi function q(x) which differs only weakly (at T close to 
T c ) from its maximum value q(l). 

There is another family of spin glasses that are characterized by one-step 
replica symmetry breaking (Parisi function in this models is a step func- 
tion [17, 10, If]). Moreover, these models are known to possess dynamic 
(at T = Td) and static T = T c phase transitions at different temperatures 
(Td > T c ). Therefore, one expects all effects of history-dependence and non- 
ergodicity to be more dramatic in the models of this second family. Note 
that the configurational entropy is known to be finite in these models right 
at the transition point: 

S con (T)/N<xQ(T D -T), (2) 

see [8] and references therein. Indeed, as it will be shown below in the 
present paper, the slow cooling solution for the model of that kind differ 
qualitatively from the results of static replica theory. Another reason to be 
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interested by this second family is the development of glass models without 
preexistent disorder [12, 13] whose behaviour seems to be similar the random 
spin glasses of the second family. Also similar dynamical equations were 
derived by Leutheusser [14] in the connection with structural glass problem. 

Recently an important progress was achieved in the investigation of the 
dynamics of spherical p-sp'm interaction model [15], which is exactly solvable 
in some sense: dynamical equations for this model can be written exactly. In 
this paper we discuss more complicated (and more realistic) case of the Ising 
p-sp'm interaction spin-glass model which belongs to the family mentioned 
above. In the case of infinite p the model is equivalent to random energy 
model and exhibits one-step replica symmetry breaking [17]. For a finite 
values of the parameter p } this phase is stable in some vicinity of the transition 
point only, whereas at lower temperatures the second phase transition with 
a full continuous replica symmetry breaking takes place [10]. This property 
distinguishes the Ising version of p-sp'm model from the spherical one and it 
is interesting to inspect how it influences the dynamics. 

The paper is organized as follows. In the sections 2 and 3 we introduce 
the model and derive averaged over disorder mean field generating functional 
for the correlation and response functions. In the section 4 ergodic dynamics 
is investigated and transition line Tr>(b) (where b is the external magnetic 
field) is found. It is shown that at high enough magnetic field b > b cr phase 
transition is of continuous nature and reminds the one known to exist in the 
SK model at finite field. The dynamic critical exponent V\ characterizing the 
decay of spin-spin autocorrelation function at the dynamic transition line, 
C(t) oc is found as function of 6, and shown to vanish as b — > b cr + 0. 

At lower external fields the transition is discontinuous: nonzero long-time 
limit of the autocorrelation function C(t — > oo) appears just at T = (6), 
at zero field C(t — > oo)\t d = p — 2 at small p — 2. The "tricritical" field 
value separating the regions of first- and second-order spin glass transitions 
is identified as b cr oc p — 2. The reason for an existence of tricritical field b cr 
is rather simple: in a non-zero external field spins are polarized already at 
T > Td(6), i.e. local < <7; >^ 0; the interaction between a = cr 8 — < Oi > on 
different sites i,j will now contain usual pair-wise random term a^jJij with 
the relative strength oc 6, which tends to produce usual SK-type transition 
and competes with an original p-sp'm interaction. Thus, at high enough fields 
the "induced" interaction wins and transition is of continuous nature. Phase 
diagram in the T — b — p coordinates is shown schematically in the figure 1. 
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Figure 1: Phase diagram. Thin lines: continuous transition, thick lines: 
discontinuous transition. 

In the section 5 slow cooling of the system in the spin-glass phase is con- 
sidered. It is shown that the behaviour of the nonergodic dynamic response 
is qualitatively the same in the whole low-temperature phase, i.e. there is 
no additional low-T phase transition known from the replica solution [10]. 
As the temperature goes down, the solution for the anomalous correlation 
and response functions interpolate smoothly between the one characteristic 
for spherical p > 2-spin model [15] and the one for the usual p = 2 Ising 
model [2]. We also found the downward jump in the specific heat as temper- 
ature decreases through T = Tjy. Such a behaviour is known to exist in real 
glasses; in our model its origin may be associated with an abrupt drop in the 
conhgurational entropy (cf. 2) due to the freezing of the system in the one 
of all possible (oc e Scom ) metastable states. 

Section 6 is devoted to the discussion of the results. 



Ising p-spin interaction spin-glass model is described by Hamiltonian [17, 10]: 



2 The model 



N 



H 





(3) 



l<t"i<— <*jp<JV 
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where b is an external field. The spin glass, described by this Hamiltonian, 
is a system of N Ising spin interacting via randomly quenched infinite range 
interactions Ji 1 ...i p . We will take, for simplicity, distribution of constants 
■h x ...i v to be Gaussian : 



P(Ji 



ii...ipj 



irJ 2 p\ 



1/2 



exp 



,Pp\ 



1 < k < ■ ■ ■ < i P < N . (4) 



We will assume Glauber dynamics for <7;: the probability for Oi to change 
it's sign during unit time T _1 is 



w 



^1, 



■ ^n)} = -(1 -<?i tanh/3/i 8 (t)) (5) 



dH 



where local field hi = — — , (3 = T 1 is inverse temperature. In the following, 

Ocr t 

we will put r = 1 without loss of generality. 

The quantities of our interest are the average response function 



G(M') 



S(3b t {t') 

that vanishes for t < t', and the average spin correlation function 



C(t,t') = (a t (t)a t (t')) 
Here (...) means the dynamic average and 777 — the average over disorder. 



3 Averaged self-consistent generating func- 
tional 

In the thermodynamic limit N — > oo, the analysis simplifies and dynamics of 
the system can be described by a set of self-consistent equations for a single 
spin. These equations can be derived by introducing a generating functional 
for Glauber dynamics [19, 20] and averaging it over disorder. This functional 

is defined as Z[X] = l^exp ^J2 J Xi(t)ai(t)dt^ ^ and, as shown in Appendix, 

can be written as: 

Z = JZ 

h=l3b 
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Here 



j = exp I J dtJ H .. Ap 8 H (t)8 t2 (t) ■ ■ ■ 8 tp (t) 



X 8 8 5 

Of. = — — , o t - 



8X(t) ' Sh(t) ' 
Z [X] is the generating functional for non-interacting spins: 



Z [\] = exp ^2 / \ % {t)m t (t)dt 



V =1 to / 
and rrii obey the equation: 

d t m t (t) = i\ t (t)(l - m1(t)) - (m t (t) - tanh/j), m t (t ) = m . (6) 

It is convenient to use this functional in another form. Note that 

Z [X] = J JlVaiVaiVhiVhiexp J X t cr t dtj exp(S) , 

i=l \i=l / 

where 

S = J - mi)dt - i J h t (h t - (3b) c/t j , 

d t m t (t) = icr t (t)(l - m 2 t (t)) - (m t (t) - tanh/36), m t (t ) = m . 
Action of operator J gives 

Z[X] = J \{Vcr % Va % Vh % Vh % exp |^ J \ % a t dt \ exp(S) (7) 



8 = 1 \i=l 

where 



S = Y J °i{ a i - rrii)dt - i j 'hi ^ - /?^rj dt j . 

This integral is normalized to unity (i.e. Z[X = 0] = 1), so that average 
over Ji 1 ...i p can be performed. The following derivation of the mean field 
generating functional is standard and almost coincides with the analogous 
one derived by Kirkpatrick and Thirumalai [11] for the case of Langevin 
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dynamics. There is no need to report this derivation and we write only 
ultimate result: 

Z[X] = I ' VhVhVoVo exp(S) (8) 



S = —i J h(h — j3b) — J [ia(a — m) — Act] + S g i , 
S gl = -J 2 ( y i i(l) l i(2)C p - 1 (l,2)h(l)h(2)+ 



+ J 2 J(p- 1)^(1)^(2)G(1, 2)C p - 2 (l, 2)»A(1)<7(2), 

dtm(t) = ia(t)(l — m 2 (tj) — (m(t) — tanh/j), m(t ) = m 



(9) 

(10) 
(11) 



with the obvious abbreviation of time arguments. Here we use the notation 
fi = p(3 2 J 2 /2. The correlation and response functions C(t,t') and G(t,t') 
have to be determined self-consistently from (8) with 



and 



C(M') 



G(M') 



S S 

J s_ 

8X t S(3b t , 



Z[X] 



Z[X] 



A=0 



A=0 



In the case of p = 2 the action in (8) coincides with the action of Som- 
mers [19] for the Sherrington-Kirkpatrick model. Further we will be mainly 
interesting in the case of e = p — 2 <C 1, since we already can consider p as 
a continuous variable. We will put also J = 1 to simplify expressions. 



4 Dynamics above the transition line 

In this section we study the dynamical mean field equation for constant tem- 
perature and magnetic field. Under these conditions, the system is expected 
to be in the paramagnetic phase (see the criterion below), so that the corre- 
lation and response functions depend only on time difference and are related 
by the fluctuation-dissipation theorem (FDT) 

G(t) = -e(t)d t c(t) . 
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Here it is supposed that t from (6) is equal to — oo, so that m does 
not depend on m . In the non-zero magnetic field, Edwards- Anderson order 
parameter q = lim^oo C(t) is not equal to zero, too, at any temperature, so 
it is convenient to represent the correlator C(t) as C(t) = C(t) -\- q. The part 
of S g i containing q is 



exp 



1 



fiq p - 1 h(l)h(2)dld2 



exp iz / h(t)dt 



where 



means 



dz 



exp 



Analogous to previous section, one can write Z[X] in the following form: 



Z 



JZ 



h=pb+z 



where J is 
J 

and 



exp 



C(l,2) 



s s 

Sh(l) 8h{2) 



+ G(l,2) 



s s 

Sh(l)8X(2) 



:i2) 



C(l,2) = ^(1)^(2)0^(1,2)- ^(1)^(2) q 
G(l,2) = {p-l)Jfi(l)fi(2) ^- 2 (1,2)G(1,2) . 



p-i 



The function m in Z obeys equation (6). 

One can show [19, 20] by means of FDT that the static limits of all 
correlation functions are independent of the short-time parts of C and G 
in (12). This leads to the equation 



m 



(13) 



where m = tanh(z + fib) and (. . .) means average over z. This equation 
coincides with the replica-symmetric equation, found in [17, 10]. 

Moreover, we will see that for t ^> 1 the expansion of C and G in (fiJ) 
from (12) is at the same time an expansion in G{u) — G'(cu)| w= o. Therefore, 
one is able to use perturbation theory with respect to C and G. 
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Following Sommers [19], we write several terms of this expansion applying 
FDT to each term: 



(1-m 2 ) iu((l-m 2 ) 2 ) 



1 — ito 



+ 



'I - ito) 



C(t)dt + C(t)\ t=0 



to 



2\3 



(1 - m 2 ) 



lid 



(1 - ILO) 



{l-lL0f 

2m 2 (l - m 2 ) 2 



C{t)dt + C{t)\ t=0 



+ 



C(t)\ dt+(C(t)\ t=0 



(14) 



Let us consider, at first, sufficiently large magnetic fields where the transi- 
tion is expected to be of the second order (the exact criterion will be derived 
below). The correlator C{u) at u — > diverges at the transition line where 

MP - 1) ((1 - m 2 ) 2 ) qr 2 = I • (15) 
At T = T c we have C(t) = At~ Vl where V\ obeys the equation: 



7T COt TTV! ((1 - m 2 f) ^4 = M 1 " m2 f) + 



(16) 



((1-m 2 ) 2 ) 



AjiqP- l {p - 1) 

For the case of p = 2, this equation was derived in [18, 19]. On the 
contrary to Sherrington-Kirkpatrick model with p = 2, v\ becomes zero at 
tricritical point where 

2 



1 



m 



,2\3 



2m 2 (l - m 2 ) 2 ) + ((1 



m 



,2\2 



P 



;i7) 



2fiqP~ 1 (p - 1) 

This equation along with equations (13) and (15) determines T tr , qt r and 
b tr . Near tricritical point we can expand left hand side of (16) over v\. This 
expansion does not contain first order term and begins with u 2 } so that v\ at 
b — b tr <^ 1 is: 

v x oc \jb - b tr ■ (18) 

We can also find how b tr tends to zero at the limit p — > 2. The expansion 
of the equations (13), (15) and (17) over m 2 yields 



b ir £ = P 



(19) 
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At the field b < b tr and T = T c the correlation function C(t) does not 
vanish in the limit t — > oo what makes phase transition to be discontinuous. 
Let us consider the case of b = 0. At u <C 1 



G(uj) = 1 + iuj e^C' p - l {t)dt + 



iuj \ e^C' p - l {t)dt 



(20) 



All results at b = will be correct only for e = p — 2 <C 1, since we restrict 
ourselves to several terms of expansion of (12). Let's assume that C(t) have 
time-persistent part: 

lim C(t) = q 

Equation (20) yields (in the main order of q) 

q - fiq 1+£ + q 2 = 

This equation has trivial solution q = and nontrivial one starting from 
some ji = ji c . This condition determine transition point: 



'/ = e, 

li c = 1 — elne + e. (21) 

These results will be found below by adiabatic cooling method. It should be 
emphasized that the transition temperature is higher than the one derived 
in the [10] with the help of replica method. 

To determine critical behaviour of function C(t) at long times let us 

/•oo 

consider it's Laplace transform C(uj) = / e lujt C (t)dt. We suppose that at 

Jo 

T slightly above T c function C{u) have a pole contribution and a remaining 
part: 

c ^ = W^lZ + e ^- < 22) 

where relaxation time r diverges when T — > T c . At the transition point 
C(uj) oc u V2 ~ l (and C(t) oc t~ V2 ). To determine z/ 2 let us consider equa- 
tion (20) at T = T c . The terms of the order of u 1 " 2 compensate each other in 
agreement with (21). The comparing of the coefficients in terms of the order 
of to 2 " 2 gives the equation for u 2 . In the main order in e 

r(l - 2v 2 ) = 2r 2 (l - i/ 2 ), v 2 w 0.395 . (23) 
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We can also find how relaxation time r depends on the temperature T 
near T c . Equation (20) at u = yields 

1 = (1 + e)nq e - 2q 

As will be shown below characteristic time f of the function C(t) is much 
smaller than r. At the interval 1/r <C u <C 1/f one can obtain in the main 
order of e: 

= /i c + — C + -C(uj ) + UujC(uj)) +— e tujt C 2 (t)dt . 



Such an equation (at e = 1) was studied by Leutheusser [14] in his description 
of the dynamics near the liquid-glass transition. Asuming that f oc (ji c — ji)~ a 
and t oc (/i c — one obtains 

1 

° " ^ (24) 
^ 2z/ 2 

It should be mentioned that for the spherical p-spin interaction model [15] 
the response function is 

G(uj) = ( 1 - iuj 1 e lujt C' p - x {t)dt\ . (25) 



o 



For small e, this expression is identical to (20). Therefore, the asymptotic 
behaviour of the functions C(t) and G(t) is the same as for the spherical 
model (at least in the main order of e), as one can see at equations (21), (23) 
and (24). 



5 Slow cooling 

5.1 Adiabatic equations and transition line 

Probably the most direct way to investigate the behaviour of the spin glass 
on a finite time scale is based on slow cooling starting in the ergodic high 
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temperature phase. We will assume that temperature and, maybe, magnetic 
field vary on a time scale of the order of t ^> T _1 . The situation is sufficiently 
complicated because the correlation and response functions depend now on 
both time arguments and no longer on the time difference. It is convenient 
to divide this functions into "fast" and "slow" parts: 



The functions C t (t — t') and G t (t — t') decay on the time scale t <C t and 
represent the dynamics within one of the "pure" states in the system phase 
space. The relevant time scale of the functions q(t } t') and A(t,t') is t . It 
turns out that it is possible to obtain the closed system of equations for the 
"slow" parts of the correlation and response functions. For the first time it 
was proposed by Ioffe at all in [3, 2] where the Langevin dynamics was used. 
Alternative method which starts from Glauber dynamics was developed by 
Horner and Freixa-Pascual in [5]. It can be proved (see Appendix 2) that 
both methods give identical equations. Here we prefer to use the second one. 

Now let us explain briefly the main idea. As it was mentioned in the 
previous section, the functions q and A at t — t' ^> i are independent of the 
short time functions C t (t — t') and G t (t — t'). Thus we can replace C and 
G by q and A respectively in the functional (12). Moreover, we also can 
consider rait) to be equal to tanh(/36). The existence of the terms d t m(t) 
and ia(t)(l — m 2 (t)) in the equation (11) leads to the fact that correlation 
functions are not equal to the asymptotic values and have also relaxation 
parts. For example, the correlator C(t) = (cr(t)cr(O)) of free spins in the 
constant field is C(t) = tanh 2 (/36) + (1 - tanh 2 (/36))e"l*l. Without these 
terms, we would obtain C{t) = tanh 2 (/36). Actual form of the time-dependent 
correlation function for the interacting spins is, of course, different from 
simple exponential relaxation; however, its particular form is irrelevant for 
the derivation of the slow cooling equations for the slow functions A(t, t') and 



If we take into account this remarks, the generating functional (12) in the 
adiabatic limit can be written as 



C(t,t') = C t (t-t , ) + q(t,t') , 



G(t,t') = G t (t-t') + A(t,t') . 



q(q,q'). 




(26) 
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J = exp 



(p-l)J ^(lH2)A(l,2y- 2 (l,2) — 



^(1)^(2) 
5 5 



^(i)a(2) 

As in the previous section, we introduce the functions 



(27) 



T=T c +0 , 



g(M') = ^fm^Fy-^t^t') - fx^- 1 

A(t,t') = {p-l)y/fi(t)fi(t')^- 2 (t,t')A(t,t') 



and rewrite (26) as 



Z[X] = { J exp / X(t)tanh h(t)dt 



h=pb+z , 



J = exp 



1 



2 J ^ 1,2 '8h(l)8h(2) 



8 8 r - . . 8 8 
+ / A 1,2) 



^(i)a(2) 



f28) 



(29) 



i^ c q p 1 \t=t c +o ■ 



The self-consistence conditions are 

q(t,t') = 8 t 8' t Z[\] ,A(t,t') = 8 t 8' t Z[\] 



where 



8, 



8 



8X(t) 



8 



8(3b(t) 



To derive equations mentioned above we have to expand J (29) and leave 
several terms of this expansion. This method is correct near the transition 
point at b = and e <C 1 and also at b near or greater than b tr and arbitrary 
e. At this conditions, the values of q t ^ and / A tjt / dt' is small. The expansion 
of J to the second order in q and A results after some algebra in equations 
(30) and (31). The essential remark should be done before we write this equa- 
tions. After the action of operator / A(l,2)$i$ 2 on exp (^f X(t) tanh f3h(t)dtj 

it appears the integral / A(t,t)(l — tanh 2 (/36)) dt containing undefined value 

A(t,t). Let us consider what such terms in the exact (not adiabatic) corre- 

r 

lators correspond to. To derive corresponding terms the action of -r*- should 

8 a 
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be extended only on the function X(t) in the exponent. We obtain as a result 
/ A(t,t')exp(t — t') dtdt' . This integral equals zero for any value A(t,t'). 



t'yt 



Therefore, in the adiabatic limit we should set A(t,t') equal to zero. 
Thus, the equations for A mentioned above is 

A t)t t = (m' 2 ) A ttt > + (m" 2 ) + -A ttt >(q t + q t >) {mm") + 



+ (m' 3 ) J A M A M / d\ + (m m'm" ) A tjt , / (A M + A ) d\ . 
The equation for q is 

1 



(30) 



qt,t> = (m' 2 ) qt,t> + - {m m" ) (q t + q t >) + (m 2 m' ) / (A + A^) dl + 



+ (m" 2 ) 



1 



1 



+ l - {m'm'" ) q tit ,(q t + q t ,) + \{m m< 4 )) (q 2 + q 2 ) + 



+ ( m 



J 3 



A M &/,i + A tl)1 q t)1 ) dl + - {m"'m 2 ) / (A t)1 q t + A t , tl q t ,) dl + 



+ (m m'm" } \ - / [A M (<£< + qt) + A t > tl (q t + qt)] dl + 



+ / {A t)1 q t)1 + A tl)1 q tl)1 ) dl + / <M A *,i + At/,0 </l ^ + 
+ i (m"m 3 ) / (At,iAt, 2 + A A ) c/lo?2+ 



+ ( m' 2 m 2 



H 



ere 



A M At/, 2 dld2 + / (At,iAi, 2 + At/,iAi, 2 j 

d n 



(31) 



m 



tanh(z + /36) . 



Several remarks should be done about the behaviour of the functions 
q(t,t') and A(t,t'). If the parameter q varies smoothly on the cooling tra- 
jectory, these functions vary also smoothly in the time scales of the order of 
to- If q jumps discontinuously at the critical temperature, we should expect 
also discontinuous jump in the function q(t,t') at t' ~ t c . This means that 
the correlation function C(t,t') varies at t' ~ t c in the time scales less than 
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the characteristic cooling time t . Thus the functions q(t } t') and A(t,t') at 
t' near t c should be related by generalized FDT: 

d t ,q(t,t') = A(M') (32) 

and consequently 

d t ,q(t,t') = A(t,t') 

what is consistent with the equations (30) and (31). This condition should be 
also satisfied if t — t' <C t . The same assumptions were used in [15] where the 
dynamics of the spherical model considered. In another words, the functions 
G(t,t') and C(t,t') at t — t' <C t and t — t' <C t c can be considered like 
^-function and ^-function respectively. 

Now we show how one can obtain with the help of this generalized FDT 
and equations (30) and (31) some of the results of the previous section. Let 
us find, for example, tricritical point. Suppose that the cooling trajectory 
cross transition line near this point at b < b tr . If we put t = t' = t c + 

q = (m 2 ) + (m /2 ) q + q 2 ((1 - m 2 ) 2 (3m 2 - 1)) . (33) 

Above T c q satisfied equation q = (m 2 ) and at b = b tr jump of the order 
parameter 8q becomes zero. Expansion of (33) in 8q gives in the first order 
marginal stability condition (15) and in the second order — equation (17). 
Consider now the case of b = and t = t' = t c + 0. The equations will be 

1 = (p- l)tiq p - 2 {l - 2q) 
q = q(l - 2q) + q 2 

This yields q . ji c = 1 + e — eloge as in the above (21). 

5.2 Anomalous response function at b = 

Let us now consider the solution for q and A in the spin-glass state below 
the transition point but at zero external field b. The equations (30) and (31) 
reduce to 

A*,*/ = A M / (l-q t -q t >) + J A M A M , dl , 

qt,t' = qt,t> (1 - qt ~ qt') + qt,t c qt',t c + (A t ,iqt',i + At/,i<2t,i) dl (34) 
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Here and further integration interval does not contain t c . The solution at 
t > t', fi = ji c + 2r and r <C e <C 1 is 

q(t, t') = (e + T + T')0(t - t'), A(t, t') = (e + r)£(t' - t c ) (35) 

At r, r' > e 

g (t, t') = VrT'0{t - t'), A(t, f) = y/re8(t - t c ), q(t, t c ) = ^/re . (36) 

As one can see, anomalous response function A consists of the ^-function 
contribution only. Such a solution, as will be discussed below, is in agreement 
with the one-step replica simmetry breaking solution obtained previously for 
the same model within the static approach [17, 10]. However, in turns out 
that these results holds only approximately and more accurate calculations 
lead to an appearance of the regular (smooth) part in the anomalous response 
function. 

To find smooth part of anomalous response function A we should expand 
J in (29) up to third order in q and A. We also can put A oc 8{t — t c ) in the 
third order terms. Thus, equation for A: 

At,t' = A t ,f [l - qt ~ qt> + 2(% 2 + ql + qlt>) + Mf- 

q-L ~ &,tc ~ 2 M',t c ] + / A M A M , dl . (37) 

Equation for q: 



Qt,t' = qt,t> 



1 - qt ~ qt' + 2(% 2 + q 2 t , + ^q^) + q t q t ' 



+ qt,t c qt,t c — 

{Qt + i 2 ')(Cc + $,t c ) _ Qt,tAt,t c qt,t> + J + qt',iA t)1 ) dl . (38) 

Consider region r < £ and introduce notation q(t, t') = e + r + r' '-\-w(t, t'). 
Then 

eA(t,t')w(t,t') + J A M A M / dl = , 

^w 2 (t, f) + J KiA + w tljl A tjl ) dl + yrr' = . (39) 

The solution of this equation has a degeneracy: we can change sign of func- 
tions A and w simultaneously. However, we should choose positive value of 
A since it corresponds to higher magnetization: 

a(m') = (6+T)S(t'-t c )+jyr £ ^-e {t -t>), w (t,t') = -H| r _ r '|^. (40) 



17 



The above choice of solution is based on the following physical argu- 
ments: the very presence of the anomalous response function is [2] due to 
the possibility for the system to choose between different configuration-space 
valleys (which come into existence during the cooling) in order to minimize 
the free energy of the final state. In particular, when cooling procedure is 
done at some constant magnetic field, the valleys with higher magnetization 
along the applied field direction are certainly preferable — which means that 
irreversible contribution to the susceptibility should always be positive. 

The result (40) can be compared with the analogous one in the spherical 
model [15]. There the smooth part of the function A at small r and e is 
also a constant proportional to yje. However, there is the family of special 
cooling trajectories where A have no smooth part and cooling at zero field 
belongs to this family. It is not clear whether such a family exists in the Ising 
model. 

Consider now anomalous response function at r, r' ^> e. The equations 
for q and A in this region are: 



A(t, f) [e log + ql t , - r 2 - r' 2 j + J A M A M , dl = 

%t> I £ log —j= + -q tit , -t -t I + 



qt,t c qt',t c + J (<?t,iAf,i + q t ',iA ttl ) dl = , 

elog-%) + / A M9Mc dl = . (41) 



qt,t c 

At r, t 1 ^> yje the terms containing e in (41) are small and solution is 
closed to one for the Sherrington-Kirkpatrick model: 

A(M') = 2r'^0(f-f'), g (M') = 2r . (42) 
In the inverse limit r, r' <^ yje we can introduce the notations A(t,t') = 

lf£ / TT 1 

J — ®(t,t') and q(t,t') = y w(t,t') and write these equations as 

0(M>(M') + / ©t,i©M' dl = , 
^w 2 (t, f) + J KiOtM + uv,i0 M ) dl + 2 -tt' -t 2 -t' 2 = . (43) 
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We have not found analytical solutions of the equations (43), but it seems 
evident that these solutions (which may be found numerically) interpolate 
smoothly between the results (40) and (42) corresponding to the regions 
(r, t' <C e) and (r, r' ^> \/e) respectively. In another terms, there is no 
qualitative difference between the solutions in these two regions. This re- 
sult looks quite surprising: it was found within static replica approach by 
Gardner [10] that the second phase transition should exist in this model (at 
t ~ \/e) } which is characterized by full replica symmetry breaking. In the 
dynamic approach we have not found room for such a transition. 

5.3 Susceptibility and heat capacity 

Now it is easy to find observable values such as field cooled and zero field 
cooled susceptibilities and heat capacity. Let's consider the region r < e. 
The value of \fc is 

Xfc= f G(t,t')f3(t')dt' = ^ 

J —CO 1 c 

with the accuracy of Ty/e and hence have a jump in the derivative with 
respect to temperature. Zero field cooled susceptibility Xzfc is determined 
by the integral of "fast" part of response function: 

Xzfc = f (G(t, t') - A(f , t')) (3(t')dt' = l ~ £ ~ T 

•J — oo 1 c 

and has a jump at the transition point. 

Heat capacity C = dE/dT also has a jump, in order to derive it the inter- 
nal energy should be written in terms of functions C(t,t') and G(t,t') [15]: 

fj(t)E(t) = - f C p - l {t,t')G{t,t')dt' 

Derivative of this equation with respect to temperature yields a downward 
jump in the heat capacity 

AC = Ct c -o — Ct c +o = —£ ■ 

(note that in the standard Landau theory of second-order phase transtions 
heat capacity shows an upward jump on temperature decrease). Static the- 
ory [11] also predicts a jump in heat capacity, but at lower (static transition) 
temperature. 
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6 Conclusions 



We have studied dynamics of the p = 2 + e-spin interaction Ising spin glass 
above and below the dynamics transition temperature (implicitly defined 
by Eq.(21). The discontinuous transition is found at zero and weak external 
magnetic fields, b < b tr ~ e, whereas at higher fields b > b tr the transition is 
continuous and resembles (though definitely is not identical to) the SK model 
transition. Near the "tricritical point" b = b cr the dynamic exponent deter- 
mining the rate of long-time relaxation right at the transition is approaching 
zero: v\ oc \Jb — b tr . 

In the glassy phase the history dependence is described quantitatively by 
the anomalous response and correlation functions A(t,t') and q(t,t'). We 
have derived equations for these functions and (in the case of zero exter- 
nal field) solved them in several regions of reduced temperature r <C 1. Very 
close to Td, at r <C e, the main contribution to the anomalous response func- 
tion A(t,t') comes from the ^-function term (cf.Eq.(35)), in agreement with 
replica theory solution [10]. Indeed, the one-step replica symmetry breaking 
found in [10] is commonly interpreted [8] mature of an instantaneous 

"appearance" of an exponentially large number of valleys right at the tran- 
sition (i.e. of the extensive configurational entropy S con oc N), those states 
do not acquire additional "fine structure" (and their number does not grow) 
as temperature decreases further. Such a structure of equilibrium valleys 
would precisely agree with the "c/e/ta-function only" solution for A(t,t'): the 
physical idea behind the slow cooling approach is that anomalous response is 
nonzero when the number of valleys grows (with T decrease usually) making 
it possible to lower free energy by the proper choice of the valley. 

However we found that even in the smallest-r range the A(t,t') response 
contains also the smooth part (cf.Eq.(40)) corresponding, within the same 
logic, to the continuous splitting of valleys below Tp. The same conclusion 
was reached by a different route in a recent preprint [8]. The above results 
provide additional (cf. [15]) evidence that the structure of valleys most rel- 
evant for slow dynamics is different from those responsible for equilibrium 
Gibbs partition function (which is reflected already in the fact that the dy- 
namic transition temperature is higher than the static one); for the discussion 
of the additional aspects of relation between dynamics and static quantities 
see [22]. 

Heat capacity and zero field cooled susceptibility have a downward jump 
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at the transition temperature. It should be mentioned that similar jump in 
heat capacity has been observed experimentally in real liquid-glass transi- 
tions [23]. 

At lower temperatures r ^> yfe the solution of slow cooling equations 
is approaching the one of the SK model [2, 4], whereas in the intermediate 
region e « t <C y/e it interpolates between the above limiting cases. Since 
qualitative features of the solutions are similar in the regions r <C e and 
t ^> y/e, we do not expect any additional phase transitions within slow 
cooling approach. Thus we found that the low-temperature properties of the 
SK (p = 2) and p = 2 + e Ising glass models are rather similar, in spite 
of the drastic difference known to exist between the corresponding phase 
transitions. 

Real experiments on as well as Monte-Carlo simulations of glassy systems 
are always done within limited time scale, which makes virtually impossible 
to observe Gibbs equilibrium properties which are the subject of static replica 
theory; on the other hand, slow cooling approach seems to be most suited to 
describe finite-time experiments. Unfortunately it does not seem possible to 
compare directly the present analytical results with Monte-Carlo simulation 
since our calculations are done for small e Cl; however we expect qualitative 
features of our solution to survive for, e.g. p = 3 Ising glass model which 
could be simulated directly. 

Let us note finally that the potentially interested problem which we have 
not studied here is the anomalous response behavior close to the tricritical 
point b = b cr (where dynamic exponents v\ and z/ 2 tend to zero). 
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Appendix 1 



Let us derive at first the generating functional Z[X] for one Ising spin in the 
arbitrary field. Distribution function for (5) obeys equation 

d t P{(T,t} + LP{<T,t} = , (44) 

where 

LP{a } t} = -J2 [w{<t <t'}P{<t, t} - w{a' a}P{a', t}] 
{a'} 

We will be also interesting in the function P(a,t,a',t') — probability to 
find spin a at the time t provided that at the time t' it's value was <r', 
P(cr, t, cr', t') = 0, t < t'. This probability is a Green's function of the operator 
d t + L a : 

(d t + La) P(<r, t, a 1 , f) = 8(t - *')<W' 
Solution of this equation is 

P(<7, t, a', f) = ^(1 + am t {a', t'j) (t > f) , 

where 

t' 

We also can write the solution of (44), i.e distribution function: 
P(a,t) = ^(1 + am t (m ,t )) 

The process described by the equation (5) is markovian. For this reason 
correlation function = (cr tl cr t2 • • • cr t „} is determined only by distribution and 
probability functions: 

C = CTlP{^l,t 1 ,(J2,t2)(J2P{(J2,t2,(J 3 ,t 3 ) X 

{a 1 a 2 ---u n } 

^n-l^K-l, tn-l, &n, t n )^ n P(^n, t n ) (45) 

for (ti > t 2 > • • • > t n > t ). 
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Later it will be convenient to introduce auxiliary functions 
A(t) = exp / A(r)m(r) dr 



and 

B(t) = m(t)exp ( j A(r)m(r) dr 

where mit) obeys equation 

d t m(t) = i\(t)(l - m 2 (t)) - (m(t) - tanh (3b), m(t ) 
This functions are connected by the relation 



B(t)= I exp (t - t) (A r + tanh /36 r )A(r) dr 
Suppose that 



(cr 2 • • • cr n ) = S 2 ■ ■ ■ S n Z 



A=o ' ^ 2 > ' ' ' > ^ n ' 



where 

Z [\] = l^exp \(t)<r(t)dtj ^ = exp \(t)m(t)dt 

and prove the same equation for n-point correlator. 
On the one hand one can find using (45) 

(a 1 a 2 ■ ■ ■ cr n ) = (a 3 ■ ■ ■ a n ) exp(t 2 - h) + 
ti 

+ (cr 2 • • • (T n ) J exp (t — ti) tanh (3b T dr . 
t 2 

On the other hand, for t > t\ > ■ ■ ■ > t n > t 



SiS 2 ---S n Z [\]\^ =o = S 1 S 2 ---S n A(t) 



A=o 



23 



S 2 ---S n B(t 1 ) 



A=o 



4 • • • 4 



exp(t 2 -t 1 )A(t 2 ) + 



exp(r — ti)(X T + tanh (3b T ) S 2 A(t) dr 



A=o 



5* ■ ■ ■ S r . 



exp(t 2 — ti)A(t) + J exp(r — t\) tanh j3b T S 2 A(t) dr 
t 2 



A=o 



which is identical to (48). 

The procedure of the generalization of the functional Z to the case of 
interacting spins is described in [19, 20]. Let us denote C = ((T\(t 2 . . . (Tk} } 
where indices are the unions of the time and space arguments. Now one can 
introduce auxiliary fields h and h and expand term with interaction over J. 
Correlator C is 



oo „ N 



N 



C =J2 / II VhiVhiexp [^2 -i [ hi(hi - (3b)dt\ (^Kn) 



n=0 i=l 



\i=l 



(49) 



where 



a x a 2 . . .a k 



52 J dt.J H ... lp ih %1 (T l2 ■ ■ ■ o lp 

\«l<«2<-"<*'p / 



and (...) means average over dynamics of the noninteracting spins in the 
field hi. (K n ) can be written variation of Z : 



52 {—,Kr}) =—,52 ■■■hi 52 \ dtj %1 ... %p ih % 2s %2 ■ ■ ■ s lp ] z 



n=0 n=0 \«i<«2< - "<*'p 



A=o 



S X S 2 ...S k exp I 52 / dt.J H ... lp ih H S l2 ■ ■ ■ S lp ] Z 

\«1 <«2<-"<«'p 



A=o 



where S t = S/SXit). In the expression (49) the integration over h and h can 
be performed, which replaces h by fib and ih by S/Sfib. Thus the generating 
functional Z is 

Z = JZ 

h=l3b 
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where J is 

j = exp I J dtJ H ... tp 8 H 8 t2 ■ ■ ■ 5i p J 

\«l<«2<-"<*'p / 

Appendix 2 

Let's show how one can derive adiabatic generating functional (26) starting 
from Langevin dynamics of soft spins: 

To^Mt) = + Wi ~ uaf + £i(t) . (50) 

Here H is the Hamiltonian (3), £i(t) is a white noise with zero mean and 
variance 

<UtMt') >=2To 1 8 ij 8(t-t f ) . 

The constants r and u should tend to infinity in the Ising limit with r /u = 1. 

As was shown in [11, 15], average this Langevin equation over disorder 
yields 

t 

_3/,\ , / j, i A/. ./\ _/ ./\ 



d t <r(t) = r (j(t) - u(T%t) + J df G(t, t')(j{f) + h(t) + £(t) (51) 

— oo 

with nonlocal noise 

< mat') >= 2T^8(t - t 1 ) + C(t, t 1 ) . (52) 

The next step is to divide the functions C and G and noise £ into slow 
and fast parts [18, 2]: 

C(t, f) = C(t - f) + q(t, t'), G(t, f) = G(t - f) + A(t, f). 

t(t) = i(t) + z(t), 

< z(t)z(t') >= y J f i(t) f i(t')q p - 1 (t,t') (53) 

Integration over fast noise leads to the equation for the slow magnetiza- 
tion [2] which replaces Langevin equation in the adiabatic limit: 

(v) z =m(t)=tznh(H eJJ (t)) (54) 
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where in the case of Ising spins 

H eJJ (t) = z(t) + (3b(t) + / dt' A(t,t')m(t') . (55) 



Now we can easily construct adiabatic generating functional Z[X\. Note 
that correlation function can be written as: 

C(ti, . . . ,t n ) = (m(ti) . . .m(t n )) z . 

Thus, 

Z[X] = Q VhVh J exp(S)exp( J X tanh (3b)^ (56) 

where 

S = -if h t (h t -z t -fib t -j A ttt , tanh h t ,dt')dt , 

and Jacobian 

d(h t -zt-pbt-f A t)t , tanh h t >dt') 



dh t , 

exp (— j A tjt (l — tanh 2 h t )dt\ = 1 



J 



since we should put A(t,t) = in the adiabatic limit, as was mentioned in 
the previous section. Average over z yields 

S = —i J h(h — (3b) — - J hih 2 yj Hi^ql^ 1 + i J h t tanh h 2 A lt2 ■ 

Then we should expand Z[X] over q and A and rewrite each term of the 
expansion as a variation of Z = exp(J X tanh h). After some algebra one can 
derive formula (26). 
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